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Abstract 

This document aims to answer the question of what is the mini¬ 
mum delay value that guarantees convergence to consensus for a group 
of second order agents operating under different protocols, provided 
that the communication topology is connected but unknown. That is, 
for all the possible communication topologies, which value of the delay 
guarantees stability? To answer this question we revisit the concept 
of most exigent eigenvalue, applying it to two different consensus pro¬ 
tocols for agents driven by second order dynamics. We show how the 
delay margin depends on the structure of the consensus protocol and 
the communication topology, and arrive to a boundary that guarantees 
consensus for any connected communication topology. The switching 
topologies case is also studied. It is shown that for one protocol the 
stability of the individual topologies is sufficient to guarantee consensus 
in the switching case, whereas for the other one it is not. 


1 Introduction 

The coordinated execution of a task by a team of multiple dynamic agents 
almost always require that a decision upon the value of a variable is reached. 
If we talk about formation flight of multiple unmanned aerial vehicles, for 
example, a common heading and a common speed must be selected. When 
this common value is to be reached by local interactions among the agents 
of the group a consensus problem appears. The work of Vicsek et al 
was the first one to study the properties of such agreement problem, con¬ 
sidering discrete agents governed by first order dynamics. This work was 
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later expanded by Jadbabaie and others. It was the work of Olfati-Saber 
and Murray the one to introduce the expression consensus. They exhaus¬ 
tively studied the conditions under which a control protocol leads a group of 
agents governed by first order continuous dynamics to reach an agreement, 
including cases with switching topologies and time delays. Many deriva¬ 
tives of this work have appeared in the years since its publication, including 
the extension to second order dynamics with the inclusion of switching 


topologies 1^, time delays |7f[14|, or both |15f[18|. A recent review paper 19 


presents a nice picture of the current challenges and results in this topic. 

One of the most important contributions of was to pinpoint the im¬ 
portance that some results from algebraic graph theory have in the topic of 
consensus. They showed that for the agents to reach an agreement upon the 
average of the initial conditions, the graph that describes the communication 
topology must be connected and balanced. Another important result is that 
the time constant of the system, i.e., the time the agent need to reach con¬ 
sensus, is related to the second smallest eigenvalue of the Laplacian matrix 
of the graph, also known as the algebraic connectivity or the Fiedler eigen¬ 
value 20 . They also show that when a time delay is present, the largest 


eigenvalue of the Laplacian limits the tolerance of the system to a delay. 
These results led to some works in which the weights of the communication 
channels are designed such that a certain delay margin is attained 

Almost all works dealing with time delayed consensus problems, whether 
they consider first or second order dynamic models, use control laws in 
which the delay affects the information of the own state of each agent as 
well as the state information coming from other agents. This so called self 
delay case can be seen, for example, in (3l[7[jl4|. The consensus protocols 
which share this structure are always described using the Laplacian matrix. 
Reference presents a different case, in which the state of the agent is 
not affected by the delay, and therefore the Laplacian matrix plays no role. 

In it is shown, however, that consensus protocols with or without 
self delay have a common structure which can be used to simplify its analy¬ 
sis with respect to the time delay. That paper shows how the characteristic 
equation of a consensus problem can be decomposed in factors of low order 
(equal to the order of the dynamics of the individual agents) with a similar 
structure. The factors differ only in a certain coefficient given by an eigen¬ 
value of a matrix related to the communication structure. The work 24 


takes advantage of this property to perform an exhaustive analysis of a pro¬ 
tocol without self delay for a fixed communication topology. In that paper 
the concept of most exigent eigenvalue.^ defined as the eigenvalue of the ma¬ 
trix which creates the most restrictive stability boundary with respect to 
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the time delay, is introduced. 

In the present paper, we aim to perform a similar analysis for a consensus 
protocol with self delay. We show how the most exigent eigenvalue in that 
case is the largest eigenvalue of the Laplacian matrix. We compare this to 
our previous result for the case without self delay and show how it can be 
used to guarantee consensus under an unknown communication topology 
under the presence of time delays in both cases. Furthermore, we study the 
switching topologies case and show how the protocol with self delay is stable 
regardless of the switching scheme, whereas the protocol with self delay can 
be destabilized if the communication topologies switches in a certain pattern. 

The paper is organized as follows. Section presents the two consensus 
protocols under study and some brief results regarding their stability analy¬ 
sis. Section presents the results on the concept of most exigent eigenvalue 
for both protocols and Section [^studies the switching topologies case. Some 
conclusions and directions for future work are presented in Section In the 
rest of the paper, scalars are represented by lowercase italic letters (fc. A, r), 
vectors by lowercase bold letters (x) and matrices by uppercase boldface 
letters (A, A). 

2 Consensus Protocols and Stability Analysis 

In this work we consider a set of n one dimensional agents driven by second 
order dynamics of the form Xi{t) = i = 1, 2 ..., n. The control 

input Ui{t) is computed based on the state of agent number i as well as 
the state of some of its peers, known as the informers of i. The set of 
informers of agent i is denoted as A/^, and its cardinality, i.e., the number 
of informers of agent i, is denoted as 6i. For this study we consider that 
the communication is bi-directional, i.e., i ^ J\fj ^ j ^ Mi. This means 
that the network can be represented by an undirected graph. The adjacency 
matrix of this graph, denoted as Ap = [a^j], is defined such that aik = 
aki = 1 when agents i and j share a communication link and aik — a^i = 0 
otherwise. The diagonal elements of Ap are considered to be zero. The 
degree matrix of the communication topology is a diagonal matrix A with 
its zth diagonal element equal to the degree of agent i. The Laplacian matrix 
of the communication topology is defined as L = A — A, and the weighted 
adjacency matrix is C = A“^Ap. 

We also assume that the communication among agents is corrupted by 
a time delay r, which is constant and uniform across the network. 

In this work we study two different control actions that can be used 
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by the agents in order to reach consensus. The first one was introduced 
in and is given as: 


Ui (t) =ki 


+ ^2 


E 


Xj {t - t) 


j&J^i 


Si 


Xj {t - t) 


Si 


- Xi (t) 


- Xi (t) 


whereas the second one, studied in 0 , is: 

Ui (t) =ki 


+ ^2 


(xj (t-r)- Xi {t - r)) 

jeAfi 

i't-x)- Xi {t - r)) 

jeAfi 


( 1 ) 


( 2 ) 


In equations Q and Q ki and k 2 are positive control gains selected by 
the user. For r = 0 both protocols guarantee consensus provided that the 
communication topology is connected. Connectivity is therefore assumed as 
granted for the rest of the paper. 

The main difference between protocols 0 and 0 lies in the presence of 
the so called self delay. In 0> agent i uses the delayed information coming 
from its informers and compares it to its own current state. In Q the delay 
is present in the state of agent i as well as in the state of the informers. This 
structural difference has some implications in the dynamics of the group, a 
fact that is presented in the following paragraphs. 

Despite the structural differences, it was shown in that the char¬ 
acteristic equations of both protocols can be expressed as a product of n 
second order factors. This factorization property stems from the fact that 
both protocols can be represented in state space as: 


X (t) = (I^ (g) Fi) X (t) + (M (g) F 2 ) (3) 


with X = [xi X 1 X 2 X 2 • * * Xn In being the identity matrix of order n and 
(g) the Kronecker product. The matrices Fi^2 depend on the particular pro¬ 
tocol and the control gains, whereas the matrix M depends on the protocol 
and the communication topology. For 0 we have: 


' 1 = 


0 1 
-ki -k2 
M = A"^Ar = C 


0 0 
ki k2 


(4) 
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whereas for (§: 

0 0 

ki k2 

M = Ar - A = -L 



(5) 


Taking advantage of the fact that the M matrix is in both cases diagonaliz- 
able, the characteristic equation of the protocols is expressed as a product 
of n second order factors: 


CE{s,t) = Ylqi {s, ki, /C2, r, A*) = 

1=1 

n 

J]det(sl 2 -Fi -AiF 2 e-"^) 


( 6 ) 


1=1 


in which A^, i = 1, 2 ..., n, are the eigenvalues of the matrix M corre¬ 
sponding to each protocol. Each one of the factors in Q corresponds to the 
dynamics of a certain linear combination of the positions of the agents in 
the group. 

The values of \i are always real for each protocol, since the weighted 
adjacency matrix C is a symmetrizable matrix and the laplacian L is 
symmetric. Furthermore, since C is a hollow stochastic matrix [^ , Gersh- 
gorin’s circle theorem guarantees that |A^| < 1 for protocol Q; whereas 
Az < 0 for protocol (§, because the eigenvalues of the laplacian are always 
nonnegative 
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Since we are assuming that the topology is connected, a special eigen¬ 
value is present in each protocol. This eigenvalue corresponds to an eigen¬ 
vector of M in which all the elements are equal. Without loss of generality, 
we denote this eigenvalue as Ai. For 0 the special eigenvalue is Ai = 1, 
whereas for protocol 0 it is Ai = 0. The factor of Q generated by Ai 
dictates the dynamics of the group decision value. That is, if all the other 
factors are stable the agents reach consensus and move together in a trajec¬ 
tory dictated by the dynamics of this factor. We call it the centroid factor. 
The other n — 1 factors define whether the agents reach an agreement or 
not, and the are called the disagreement factors. 

Reference 24 showed that for protocol © the group decision value is 
dictated by a weighted average, defined as: 


{t) = 


E •=! 


(7) 
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for which the dynamics is given by: 


+ (^25 + fci) (1 — — 0 


( 8 ) 


Notice that Q has a root at 5 = 0 for any value of r. This means that this 
factor is at best marginally stable, which implies that consensus, if reached, 
will be at a constant position with zero velocity. 


showed that the agents reach an average consensus given by: 


For protocol H, on the other hand, the previous works m and 30 


Cl (t) 


n 


(9) 


with dynamics: 


5 ^ = 0 


( 10 ) 


Equation (10) has two roots at the origin, which do not depend on the delay. 
This indicates that protocol (§, when it is stable, guides the agents towards 
a constant velocity with linearly increasing position. 

For the agents to reach consensus using protocol 0 n — 1 disagreement 
factors of the form 


Qi {s, ki, k 2 , T, Xi) = + {k 2 S + ki) (l - AjC = 0 (11) 

with i = 2, 3, ..., n, must be stable. When protocol 0 is used, the n — 1 
disagreement factors have the forn|^ 

qi (s, ki, k2, T, Xi) = + Ai {k2S - ki) e~'^^ = 0 (12) 


Factors 0 and ( [T^ are stable for r = 0. As r increases, the character¬ 
istic roots of the factors move in the complex plane until the delay reaches a 
critical value for which a root crosses to the right of the complex plane, intro¬ 
ducing instability. The Cluster Treatment of Charaeteristie Roots Paradigm 
(CTCR) is used to detect these destabilizing crossings. In order to limit 
the length of this paper, the details of the deployment of CTCR to the fac¬ 
tors of the form 0 and 0 are left out. They were already presented 
in (for 0) and (for d^). We only state here the final results. 

For factors of the form (|11[), the first destabilizing root crosses at a 
frequency: 

= ki- ^kl + ^ (13) 


^Notice that a change in the sign of Xi has been introduced in (12). This is to use the 
eigenvalues of L, which are all positive, and not those of — L. 
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with: 


/i = 1 — Af, and 7 = — 2fci)^ — 4fc^/i (14) 

and this crossing occurs at a delay value which depends on the sign of Xf. 


r = - (^arctan j j for A, > 0 

r = yarctan( ^,_^ ~7_p) )+") f" ^< < » 

For protocol (§, the expressions for the crossing frequency and delay 

^ kiX^ + JklXf + AkfX] 

^ ~ 2 

1 / k2UJ 

T = — arctan —— 

(jj \ ki 


(16) 


are: 


(17) 

(18) 


The stability analysis technique presented in departs from the knowl¬ 
edge of the communication topology. From here, the eigenvalues of the 
weighted adjacency or the the laplacian matrix are obtained, and the stabil¬ 
ity region in the parametric domain (/ci, /c 2 , r) is obtained for each factor. 
These regions are then intersected to find the combinations of parameters 
which make the complete system stable. 

In the next section, we show that an a priori knowledge of the commu¬ 
nication topology is not needed to find a delay bound. 


3 The most Exigent Eigenvalue 

Section [^showed that each eigenvalue of M introduces a stability boundary. 
However, there is always one eigenvalue that introduces the most restrictive 
of these boundaries with respect to the time delay, therefore defining the 
global stability region. This eigenvalue is declared the most exigent eigen¬ 
value. This section formally defines this concept, first introduced in p4| , 
and declares which eigenvalue is the most exigent for protocols 0 and^^. 

Definition 1. For a group of agents interaeting under any of the eonsensus 
protoeols defined in Q or 0. the most exigent eigenvalue is the eigen¬ 
value of the eorresponding M matrix that generates the partieular faetor 
Qi (s, /ci, ^ 2 , T, Xi) in the eharaeteristie equation 0 whieh introduees the 
smallest destabilizing erossing as r inereases starting from 0, for a fixed set 
of ki and k 2 values. 
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Lemma 1. For the consensus protocol without self-delay defined by 0; the 
most exigent eigenvalue is the smallest eigenvalue of the weighted adjacency 
matrix C. 

Proof: The proof was presented by . □ 

Lemma 2. For the consensus protocol with self-delay defined by 0; the 
most exigent eigenvalue is the largest eigenvalue of the laplacian matrix L. 

Proof: We prove this following the same path used by to prove 

Lemma [Tj We show that the following statements are true when we consider 
cj > 0: 


(a) (jj is a, monotonic one to one function of A and dujjdX > 0. 

(b) T is a monotonic function of uj and dr jduo > 0 . 

(c) T is a monotonic function of A and dr/dX < 0. 

In order to prove (a), we use = 7 and consider 0 for two different 

A, recasting the equations as: 


7 " - (A: 2 Ai )7 - (^lAi)" = 0 (19a) 

- {k2\2? 1 - iki\2? = 0 (19b) 


If this two polynomial equations have a common root, the Silvester resultant 


matrix: 


R 


-{k2\ir 

1 

- {k2\2? 

1 


(kiXif 

(/C2 Ai)^ 

{kiX2f 

{k2X2f 


0 

-{k^X^f 

0 

- (A:iA 2 )" 


( 20 ) 


must be singular. We have that det (M) = “ ^ 2)^5 which can be zero 

if and only if Ai = A 2 . Therefore, 7 is a one to one function of A. Since 
-^ 7 , the first part of (a) is proven. 


UJ 


By using implicit differentiation, we see from (19a) that: 




(W 

27 


( 21 ) 


and combining ( 21 ) with du/d'y = 1 /( 2 .^) we have: 


duj duo dy {k 2 X)^ {k 2 X)^ 

dX dy dX 47 ^/^ 

8 


( 22 ) 









which proves the second part of (a). 


To prove part (b), we consider (18) and take its derivative with respect 
to oj: 

dr k 2 /ki 1 f k 2 UJ \ 

- = - — -I - I 


duj 


CO 1 + U¥^ 


( k20j\ 

I fcl ) 


LO^ 


\h J 


(23) 


We are interested in showing that dr jduj < 0. By multiplying (23) times 
we obtain the following inequality: 


k2(jO 

ki 


1 + 


V ) 


f k2(jj\ 

7 T — arctan —— < 0 

^ V / 


(24) 


Using the substitution = tan(0), valid for 0 < 0 < 7r/2, we transform 
([^ into: 


tan0 


1 + tan^ 9 


9 <0 


(25) 


which can be simplified as: 


- sin (29) — 9 < 0 


(26) 


Figure [^presents a plot of the left hand member of (26), which clearly shows 
that the inequality is true in the interval of interest. This proves part (b): 
dr jduj < 0. 

Finally, statement (c) is a direct conclusion of (a) and (b). Since r is a 
monotonic function of uj and uj a monotonic function of A, using the chain 
rule we observe that: 

dr dr duo , 


which implies that the largest value of A always invites the minimum r, and 
therefore the most restrictive stability boundary. This completes the proof 
of lemma [2l □ 

The knowledge of which eigenvalue creates the most restrictive value 
implies that one does not need to know under which communication topology 
the agents are operating in order to guarantee stability with respect to the 
delay. 

Let us consider a group of agents operating under protocol 0. Accord¬ 
ing to Lemma [T| the most exigent eigenvalue is the smallest eigenvalue of 
the weighted adjacency matrix of the communication topology under which 
the agents are operating. But this eigenvalue is lower bounded: it can never 


9 










Figure 1: Plot of f{9) = ^ sin(20) — 9 for 0 < 9 < 7r/2. 


be smaller than —1. This implies that the stability boundary defined by 
([l^ for A = —1 is the most restrictive for any possible topology. Figure 
^ shows a plot of this boundary in the tree dimensional domain (fci, k 2 ^r). 
Selecting a parametric combination below this surface guarantees that the 
agents reach consensus regardless of the topology under whieh they are oper¬ 
ating^ provided that connectivity is present. This is also independent from 
the number of agents. 

In order to perform a similar analysis for protocol Q , we need to consider 
an upperbound for the eigenvalues of the Laplacian of a graph. The first and 
most conservative of such bounds was presented by [^, and is states that 
the eigenvalues of the Laplacian of a connected graph are always less than 
or equal to twice the largest vertex degree of the graph, i.e., A < 2max^ 5i. 
Some other authors have presented less conservative bounds imi, but they 
requiere to have an extra knowledge of the structure of the graph. Without 
any a priori knowledge of the structure of the communication topology, 
using the bound given by is the safest approach. 

For a group of n agents the highest possible degree is, of course, n — 1. 
One could then consider A = 2 (n — 1) as the upperbound for the most 
exigent eigenvalue. However, stated that the upperbound is reached 
when the graph is bipartit^ and regula]0 The highest possible degree for 


graph is bipartite if its vertices can be separated in two sets such that no vertex is 


adjacent to another member of its set 29 


29 


A graph is regular if all vertices havethe same index, i.e., Si = Sj for every i, j G [1, n] 
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Figure 2: Absolute stability boundary in (fci, /c 2 ,t) domain for a group 
of agents operating under protocol 0 with an unknown communication 
topology. Parametric selections below the surface guarantee stability and 
convergence to consensus. 


a bipartite graph with n vertices is n/2, and therefore the highest Lapla- 
cian eigenvalue is given by A = n. This is indeed supported by the results 


presented in 33,34 


Since the largest and most exigent Laplacian eigenvalue is dependent on 
the number of agents, so is the delay margin for a team operating under 
protocol Q. Figure 3(a) shows the absolute stability boundary for a case 


with 6 agents, whereas Fig. |3(b)1 shows it for a group of 10 agents. Although 
the difference is not too pronounced, it is possible to see that the boundary 
is more restricted for larger groups. From an intuitive point of view, this 
makes sense, because with more communication channels the negative effect 
of the delay is multiplied. 


4 Switching Topologies 

Section 1^ shows that we can hnd a parametric combination such that proto- 
cols 0 and 0 guarantee consensus regardless of the topology under which 
the agents are communicating, provided that it is connected. In this section, 
we focus our attention to the case in which the topologies are switching, and 
check if by stabilizing all the possible communication topologies the stability 
of the switching system is guaranteed. 
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Figure 3: Absolute stability boundary in (fci, k 2 ^r) domain for groups of 
6 (a) and 10 (b) agents operating under protocol (§ with an unknown 
communication topology. Parametric selections below the surfaces guarantee 
stability and convergence to consensus. 
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Let us consider again the state space formulation Q. In order to ob¬ 
tain a factorized characteristic equation, used a state transformation 
defined by the diagonalization of matrix M. If T“^MT = A, where A is 
a diagonal matrix, the state transformation is defined as ^ = (T“^ (g) I 2 ) x. 

The new state vector ^ can be seen as the concatenation of n state vectors 
r • iT 


dynamics: 


z = l, 2, 


.., n, each one corresponding to a subsystem with 


^i{t) 

. iiit) 



+ AiF2 


’ Ciit-r) ' 
. ii{i - t) . 


(28) 


The characteristic equation of each subsystem of the form ( |28[ ) corresponds 
to one factor qi{s^ fci, A: 2 , r, Xi) in the global characteristic equation (§. 
For i = 1 we have the system that creates the centroid factor, and it is 
therefore called the centroid system. We refer to other subsystems in (28) 
as the disagreement systems. The disagreement vector is now defined as a 
vector containing the state of all the disagreement systems. 

For a group of agents operating under a fixed communication topol¬ 
ogy, a parametric selection within the boundary defined by the most exigent 
eigenvalue guarantees that all the disagreement systems are stable, implying 
that the norms of their states decrease exponentially with time. The dis¬ 
agreement vector follows this behavior, its norm decreases until consensus 
is reached. 

Under switching communication topologies, the disagreement vector be¬ 
haves in this same way between switching instants. Whatever occurs with 
the norm of the disagreement vector at the precise moment of switching, de¬ 
fines whether consensus can be reached under switching or not. If the norm 
of the disagreement vector is continuous during all the operation, i.e., if its 
norm before and after the switching remains the same, it is exponentially 
decreasing through the operation and we can guarantee that consensus is 
reached. On the other hand, if there are jumps at the moment of switching 
and the jumps cannot be shown to be always decreasing, it is not possible 
to guarantee that consensus is reached under switching topologies. This sit¬ 
uation depends on the protocol in use, and the following paragraphs discuss 
what happens in each case. 


4.1 Protocol ([^ 

For a group of agents operating under protocol 0 , the group decision value, 
presented in is topology dependent. If the topology switches, the agents 
try to arrive to a different decision value, i.e., a jump in the norm of the 
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Figure 4: Two different connected topologies for six agents. 


vector occurs. This also causes a jump in the norm of the 

disagreement vector at the switching instant. Because of the unknown 
nature of the transformation matrix, there is no way to guarantee that the 
jumps are bounded. This behavior therefore, may induce instability, even 
when the system is switching among stable topologies. 

To better illustrate this idea, consider a group of six agents operating 
under protocol Q and switching among the two different communication 
topologies depicted in Fig. With a parametric setting of (/ci, /c 2 , r) = 
(5, 0.2, 0.06) we guarantee stability for each individual topology. We defined 
the switching pattern by means of a periodic signal with period T — 1.4 s 
and a duty cycle a, such that during a% of the period the agents commu¬ 
nicate according to the topology depicted in Fig. |4(a)[ and according to the 
topology in Fig ]4(b)| during the rest of it. 

Figurej^depicts the time evolution of the positions of the agents (Fig. |5(a)] ), 
the topology dependent weighted centroid (Fig. 5(b)[ ) and the norm of the 
disagreement vector for a case in which a = 10. It is clear that, despite the 
jumps on the value of the weighted centroid, the agents reach consensus. 

By changing only the duty cycle of the switching signal, to a = 60, 
completely different results are obtained. Figure]^ shows how the positions 
of the agents, the weighted centroid and the norm of the disagreement vector 
evolve under this conditions. 

This example shows that when a multi agent system operates under 
protocol 0 and switching topologies, it is not possible to guarantee stability 
in the general case. The conditions under which the switching system is 
stable or unstable, remain as an open question for further study. 
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Figure 5: Behavior of six agents operating under protocol 0 and switching 
stably among the two communication topologies of Figj^ a) Positions of the 
agents, b) position of the weighted centroid, c) Norm of the disagreement 
vector. 
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(a) 




(c) 


Figure 6: Behavior of six agents operating under protocol Q and switching 
unstably among the two communication topologies of Figj^ a) Positions of 
the agents, b) position of the weighted centroid, c) Norm of the disagree¬ 
ment vector. 
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4.2 Protocol Q 

For protocol (§ the situation is better. The following lemma proves the 
most important point. 

Lemma 3. For a group of agents interacting under protocol (§ the norm 
of the disagreement vector is continuous at the switching instant. 

Proof: Since the Laplacian matrix L is symmetric, the diagonaliza- 

tion T“^MT = A is performed using an orthogonal matrix T, for which 
= T^. When an orthogonal matrix multiplies a vector, the norm of the 
vector remains unchanged. Consider then a switch between two different 
topologies, described by Laplacians Li and L 2 which are diagonalized by 
orthogonal matrices Ti and T 2 . The disagreement vectors before and after 
the change in topology are given by: 

= (Tf ®l 2 )x(t) (29a) 

= (Ti’®l2)x(t) (29b) 

Given that the agents are not changing their positions or velocities when 
the communication topology changes, ||x(t)|| is the same in both equations 
( [ 2 ^ . Now, considering that (g) I 2 is an orthogonal matrix, we have that 
||^x(^)ll = ll^(^)ll = 11 ^ 2(^)115 which implies that the disagreement vector is 
continuous at the switching instants. □ 

The implication of Lemmais that if a multi-agent system operating un¬ 
der protocol Q switches among stable communication topologies, it reaches 
consensus regardless of the switching scheme. In this situation, the most ex¬ 
igent eigenvalue is very useful. It allows the designer to select a proper set 
of parameters which guarantees consensus without knowing the topologies 
or the switching scheme. 

An example of this is presented in Fig. It shows the results of a 
simulation in which six agents operate under protocol Q and switch among 
the two topologies of Fig. The switching is again periodic, and the period 
and duty cycle are the same as in the unstable example of Fig. ^ The 
parameters are selected below the corresponding surface in Fig. |3(a) Notice 
that Fig. |7(c)| shows a smooth trace, confirming the results of Lemma 

5 Conclusions 

This paper studies two consensus protocols for groups of agents driven by 
second order dynamics and affected by communication time delays. The 
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Figure 7: Behavior of six agents operating under protocol (§ and switching 
stably among the two communication topologies of Figj^ a) Positions of the 
agents, b) position of the weighted centroid, c) Norm of the disagreement 
vector. 
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agents are assumed to be operating under an bidirectional scheme, such 
that the communication topology is described by an undirected graph. The 
relation between the eigenvalues of a matrix associated to such graph and 
the maximum delay value for which the agents are able to reach consensus 
is studied in detail for both cases. 

It is shown that for one of the protocols, there is a possibility of defining 
a stability boundary that guarantee convergence to consensus even without 
knowing the number of agents. For the other protocol, the knowledge of the 
number of agents suffices to define a stability boundary which guarantees 
consensus for any connected topology under which the agents could be op¬ 
erating. The eigenvalue that defines the most restrictive boundary in each 
case is called the most exigent eigenvalue. 

The switching topologies case is also taken under consideration. It is 
shown that while for one of the protocols the stability can be guaranteed 
regardless of the switching scheme, for the other one this is not the case and 
some switching schemes may lead to instability. 

Further questions on this topic are related to another concept introduced 
in [^: the most eritieal eigenvalue. The question to be answered in this 
case is: what is the parametric combination (fci, /c 2 , r) which guarantees the 
fastest convergence to consensus? 

The codes used to create the examples of the paper are publicly available 
in https://db.tt/KeBrsfSD, or can be requested via email to the author. 
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